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, We discuss spherically symmetric static solutions of the Einstein-Klein-Gordon equations 

■ for a real scalar field with a mass and a quartic self-interaction term. As for the massless 

' case the solutions have a naked singularity at the origin. However, linear stability analysis 
shows that these solutions as well as the massless ones are dynamically unstable. 

> ' 

\^ ' 1 Introduction 

in 

' The recent developments in particle physics and cosmology suggest that scalar fields 

may have played an important role in the evolution of the early universe, for instance 
in primordial phase transitions, and that they may make up part of the dark matter. 
Q>^ I These facts motivated the study of gravitational equilibrium configurations of scalar 

' fields, in particular for massive complex fields, which form so-called boson stars 

qh| (for a review on this subject see Refs. 1 and 2). Solutions to the Einstein- Klein- 

Gordon equations for massless real scalar fields were first studied by Buchdalcl and 
more recently by WymanQ. All known static solutions for the real scalar field to 
the Einstein-Klein-Gordon equations are either topologically non-trivial or have a 
logarithmic singularity. The latter is an example of a so-called naked singularity. 
Here we present the static, spherically symmetric solutions for a real scalar field 



^ I obtained by including a mass and a quartic self-interaction term. We also discuss 

the dynamical stability of the equilibrium solutions using linear perturbation theory. 
Using the variational principle, from which the pulsation equation is derived, we 
find for all equilibrium solutions a negative upper bound for the lowest eigenvalue. 
We, therefore, conclude that all equilibrium configurations, including the massless 
ones, are dynamically unstable. 

2 Basic equations and their properties 

We consider the action of a real scalar field (j) with mass and a quartic self-interaction 
term (A > 0) minimally coupled to gravity with the following spherically symmetric 
metric 

ds^ = e^dt^ - e^dr^ - r^{d^^ + siri^Mip^) . (1) 

The static equilibrium configurations are determined by a coupled system of equa- 
tions, which can be derived from the action, by varying with respect to the various 
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fields. For all solutions we require as boundary conditions that the metric is asymp- 
totically flat and that cj) vanishes at infinity. It turns out that at the origin 4> is 
singular and behaves for all cases as: ~ —In r. For r — > 0: ^ r"" (a is real 
and a > 1) and — > 0, whereas e'' ^ [b is real and b > —1). For e'^ there are 
three different possible behaviors: e"^ = 1, e"^ ^ oo or e"^ ^ 0. See Ref. 5 for a 
plot of the results of the numerical integration for (jj and the metric functions. The 
corresponding ADM mass of the boson star has been also computed numerically. 

3 Dynamical stability 

We turn now to the problem of dynamical stability and consider small time depen- 
dent radial perturbations, which still preserve spherical symmetry. The equations 
governing the linear perturbations are obtained by expanding all functions to first 
order (A(r, i) = Ao(r) -I- 5\{r,t);v{r,t) = vq + 5v{r,t)](j}{r,t) = (j)o{r) + 5(j){r,t)) 
and by linearizing the Einstein-Klein-Gordon equations. We make also use of tlae 
Goi component of the Einstein equation, which can be integrated once in timeH. 
This way the metric-functions 5X and 5v can be eliminated from the linearized 
scalar wave equation Elu. Furthermore, we suppose a time dependence of the form 
gicrt j^Qj. §(^(^j.^iy Thus d(t){r,t) — e^'^*-'^{r)/r, where we denote the radial part of 
5(j) by \l/(r)/r. Performing a change of variable defined by ^ = g'^^o-'^o)/'^ with 
p{r = 0) = 0, the pulsation equation transforms to the following Schrodinger-type 
equation 

- + n^P])* = , (2) 

where 

V{r) = e''°-^°{{v'^->!^)/2r~A^Grcj,'S{2lT + v'^-X'^) 

-Kl67rGr</>;,e^°(m2,/)o + Xcj^l) + e^°{m'^ + SA^^)) (3) 

(r = Dynamical instability occurs whenever the lowest eigenvalue ctq is 

negative. Indeed, the perturbation 5(j) ~ e"^°* will then grow exponentially. Since 
the asymptotic behavior of the solution (t>Q{r) does not depend on the mass nor on 
the quartic self-interaction term, the behavior of V{p) will not depend on m and 
A and thus will be the same for all types of solutions. One finds V{p) ^ for 

p ^ {). It is useful to write the potential V{p) as: V{p) — — ^ + V {p) . This way we 

can write the Hamilton operator as H = Hq + V{p), where Hq — — ^ — Hq is 

not selfadjoint on its natural domain of definition 'D{Hq) — H since 

V{Hq) 7^ I?(iJg). This problem can be solved by extending in an appropriate way 
the domain of definition of Hq, such that for the correspondjng operator Ho a (a is 
real) I?(iIo,a) C V{H*) andI?(iJo,a) = V{H* J. We defineEll?(iIo,a) = 2?(i?o)+*a 
with 

= ^iH^'\expii^)p] + exp{ia)H^'\expii^)p]) , (4) 

and Hq^'' is the Hankel function. The spectrum of Ho,a depends on the value of 
a, which parameterizes the extension. For our problem V&a must be real, since it 
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describes a perturbation of the equilibrium solution for a real scalar field, 'fa is 
real only for a — 0. On this domain of definition I'(i?o,a=o) Hq (means -ffo,a=o 
from now on) has a discrete eigenvalue E = 0, whose eigenfunction v]/q=\I>^^q is in 
jC^{dp, (0,00)). From the scaling property of Hq, it follows that this eigenvalue is 
infinitely degenerated. 

Using ^'o f3{p)='^oipP) as a trial function we get upper bounds for the eigenvalue 

^2 < < ^0,/3 \Hq + V\ ■fo,0 > ^ < ^o,g I V I ^0^^ > 
~ < *0,/3 I *0,/3 > < *0,/3 I ^'0./3 > ' 

For all equilibrium solutions we easily find a value for /3 such that < \l/o,/3 | V \ 
^'o,/3 > is negative. Indeed, for large /3, \E'o(/3/c) reaches its maximum at smaller 
values of p than ioi (3 — 1. V{p) is negative at the origin and behaves as —£°i^sl for 
/9 — > 0. Thus by appropriately choosing (3, ^o(/^p) has most of its support where 
V{p) is negative, which leads then to a negative value for < ^fg.^ | V \ ^o,/3 ^■ 
Therefore, we conclude that all the static solutions are dynamically unstable 
This result is in agreement with the cosmic censorship conjectured, which excludes 
spacetimes with naked singularities. 

4 Conclusion 

From the above considerations it follows that only with complex scalar fields one 
may hope to form boson stars, and it is thus natural to study their formation and 
properties within the standard cosmological model. To that purpose an analysis of 
the coupled Einstein-Klein-Gordon equations using the Friedmann-Lemaitre metric 
has been carried out for a complex field in Ref. 9. Moreover, the time evolution 
of the perturbations of the complex scalar field and the metric has been analyzed 
in Ref. 10. It turns out, that during the oscillatory phase after inflation the per- 
turbations of the complex scalar field at best oscillate. Therefore, the formation of 
boson stars in a universe driven by the same scalar field is not possible. 
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